AC resistivity of c?-wave ceramic superconductors 
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We model d-wave ceramic superconductors with a three-dimensional lattice of randomly distributed 
TT Josephson junctions with finite self-inductance. The linear and nonlinear ac resistivity of the d- 
wave ceramic superconductors is obtained as function of temperature by solving the corresponding 
Langevin dynamical equations. We find that the linear ac resistivity remains finite at the temper- 
ature Tp where the third harmonics of resistivity has a peak. The current amplitude dependence 
of the nonlinear resistivity at the peak position is found to be a power law. These results agree 
qualitatively with experiments. We also show that the peak of the nonlinear resistivity is related to 
the onset of the paramagnetic Meissner effect which occurs at the crossover temperature Tp, which 
is above the chiral glass transition temperature Teg. 

PACS numbers: 75.40. Gb, 74.72.-h 



I. INTRODUCTION 

The interplay of superconductivity and disorder in 
granular superconductors has been of great interest, par- 
ticularly *jsffarding their magnetic properties and glassy 
behavior.S'El Granular superconductors are usually de- 
scribed as a random network of sunerconducting grains 
coupled by Josephson weak links.trQ In the last years, 
there has been a renewed interest in the study of this 
problem in high temperature ceramic superconductors 
(HTCS). Several experimental groups have found a pargu 
magnetic Meissner effect (PME) at low magnetic fields.El 
Sigrist and RiceO proposed that this effect could be 
a consequence of the intrinsic unconventional pairing 
symmetry of the HTCS of dos2-y2 typefl Depending on 
the relative orientation of the superconducting grains, 
it is possiblft_to have weak links with negative Joseph- 
son couplingaQ. These negative weak links, which are 
called TT-junctions, can give rise to the PME according to 
Ref. y,|6|. In fact, a model d-wave granular superconduc- 
tor, consisting on a network of Josephson junctions with a 
random concentration of 7r-junctions and including mag- 
netic screening, has been able to explain the paramage- 
netic Meissner effect observed at low fields E Also in this 
model, a phase transition to Syxhiral glass has been pre- 
dicted for zero magnetic fieldsJj E3 The frustration effect 
due to the random distribution of w junctions leads to a 
glass state of quenched-in "chiralities" , which are local 
loop supercurrents circulating over grains and carrying 
a half-quantum of flux. Evidence of this transition has 
been related to measurements of the nonlinear ac mag- 
netic susceptibility.li3 Moreover, the random 7r-junction 
model has also been adequate to explain several dynam- 
ical phenomena observed iiuHTCS such as the anoma- 
lous microwave absorptionjU the so called compensation 
effectlij and the effect oLapplied electric fields in the ap- 
parent critical current .E3 ■— ■ 

In recent experiments Yamao et alrJi have measured 



the ac linear resistivity, po, and the nonlinear resistivity, 
P2, of ceramic superconductor YBa2Cu408. po and p2 
are defined as the first and third coefficient of the expan- 
sion of the voltage V{t) in terms of the external current 

Iext\t)- 



V = pde. 



P2lext 



(1) 



When the sample is driven by an ac current Iext{t) — 
Iq sin(cjt), one can relate po and p2 to the first harmonics, 
V^, and third harmonics, V^^, in the following way 
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V;^^ = — I V{t)sin{jiujt)d{ujt), n = 1,3 



(2) 



The key finding of Ref. |7| is that po does not vanish 
even at and below the intcrgrain ordering temperature 
Tc2- On the other hand, p2 has a peak near this tem- 
perature, which was found to be negative. In the chi- 
ral glass phase the C/(l) gauge symmetry is not brokea. 
and the phase of the condensate remains disordered. El Ell 
The chiral glass phase, therefore, should not be super- 
conducting but exhibit an Ohmic behavior with a finite 
resistance. j_Based on these theoretical predictions Ya- 
mao et aLEZi speculated that their results give further 
support to the existence of the chiral glass phase, in ad- 
dition to previous results from magnetic susceptibility 
measurements£3 ■— ■ 

Another interesting result of Yamao et aLO is the 
power law dependence of \V^^{Tp)/ IqY\ (or of P2) at 
its maximum position Tp on Iq: \V.l,^{Tp) / Iq\ ~ ^cT"- 
The exnerimental value of the power law exponent was 

a«i.i.y 

The goal of our paper is two- fold. First, we try to 
reproduce the experimental resultslHI using the model of 
the Josephson junctions between d-wave superconducting 
grainaalj where the screening of the external field by su- 
percurrents is taken into account. Second, we discuss the 



question if the temperature Tp of the nonhnear resistiv- 
ity peak andithe transition temperature T^g to the chiral 
glass phasaly'til are related. We calculate the linear and 
nonlinear ac resistivity by a Langevin dynamics simula- 
tion. In agreement with the experimental dataO we find 
that po remains finite below and at the temperature Tp 
where p2 has a peak. Furthermore, the maximum value 
of 1^3uj(^p)/-^o I i^ found to scale with Iq with a power 
law exponent 01.= 1.1 ± 0.6, which is close to the exper- 
imental value.EII However, we find that Tp is above the 
equilibrium chiral glass transition temperature T^g- 



Here we add to the Josephson current the contribution 
of a dissipative ohmic current due to an intergrain resis- 
tance R and the Langevin noise current r/y (t) which has 
correlations 

2kT 

{ri^,m,y{t')) = -^k.'hj'&it - t') (6) 

The dynamical variable in this case is the gauge invariant 



phase difference 1 



(^i—Sj—Aij . The total flux through 



each plaquette p depends on the mesh current C„ 
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II. MODEL 



We neglect the charging effeats, of the grains and con- 
sider the following HamiltonianQJj 
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A{r)df , (3) 



The mesh currents Cp, the link currents /^ and the exter- 
nal current I^xt are related through current conservation. 
At this point, it is better to redefine notation: the site of 
each grain is at position n = [ux, ny,nz) (i.e. i = n); the 
lattice directions are /x = x, y,z; the link variables are 
between sites n and r\. + fx (i.e. link ij = link n, fj,); and 
the plaquettes p are defined by the site n and the normal 
direction /i (i.e plaquette p = plaquette n, /.t, for example 
the plaquette n, z is centered at position n -I- (x -I- y)/2). 
The current /^ (n) is related to the mesh currents C^, (n) 
and the external current in the y-direction as 



h{n) = sxf,^A C^(n) + Sx.yl, 



X.y^ext 



(8) 



where 9i is the phase of the condensate of the grain at 
the z-th site of a simple cubic lattice, A is the fluctuat- 
ing gauge potential at each link of the lattice, (po denotes 
the flux quantum, Jij denotes the Josephson coupling be- 
tween the i-th and j'-th grains, C is the self-inductance 
of a loop (an elementary plaquette), while the mutual in- 
ductance between different loops is neglected. The first 
sum is taken over all nearest-neighbor pairs and the sec- 
ond sum is taken over all elementary plaquettes on the 
lattice. Fluctuating variables to be summed over are the 
phase variables, Oi, at each site and the gauge variables, 
Aij , at each link. <I>p is the total magnetic flux threading 
through the p-th plaquette, whereas $p^' is the flux due 
to an external magnetic applied along the z-direction. 



^ext 



HS if p is on the < xy > plane 
otherwise , 



(4) 



where S denotes the area of an elementary plaquette. 
For the d-wave superconductors we assume Jij to be an 
independent random variable taking the values J or —J 
with equal probability (± J or bimodal distribution), each 
representing and tt junctions. 

In order to study the dynamical response and trans- 
port properties, we model the current flowing between 
two griSjHis with the resistively shunted junction (RSJ) 
modeljBu which gives: 
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(5) 



where ex^^i, is the Levi-Civita tensor, A is the backward 
difference operator, A~Ci,{n) = C^{n) — C,y(n — /i), and 
repeated indices are summed. It is easy to verify that 
Eq.(8) satisfies current conservation. The magnetic flux 
<i>A(n) and the gauge invariante phases 9„{n) — A+6'(n) — 
yly(n) are related as 



1>A(n) = -^e,^,A+0,(n) 



(9) 



with the forward difference operator A+6'^(n) = 6^{n + 
/i)-^.(n). 

Then, from equations (5), (6), (8) and (9) we obtain 
the following dynamical equations: 



h d0^,{n) 2e 



2eR dt 



Jfj,{n) sin6',,(n) - 5^,ylext 
^ A;[A+0^(n)-A+0.(n)] (10) 



2e£ 

-'7/.(n,t) 

which represent the RSJ dynamics of a three dimenaianal 
Josephson junction array with magnetic screening. ofl 

We can also obtain these equations from Eq. (3) if 
we add to Ti the external current term: I-Lt — Ti + 
X^n T"^e2;t^y(n). Then an equation of the Langevin form 
is obtained by taking derivatives with respect to the 
gauge invariant phase difference: 



h d0f,in) 2e SHt 



2eR dt 



.(n) 



vA^,t) , 



(11) 



leading to the RSJ dynamical equations of (10).ll3 

In what follows we will consider currents normalized 
by /,/ — 2eJ/h, time by r = (J)q/2ttIjR, voltages by RIj, 
inductance by 0o/27r/j and temperature by J/ks- We 
consider open boundary conditions for magnetic fields 
and currents in the same way as defined in Ref. W,W. 



III. RESULTS 

The system of differential equations (10) is integrated 
numerically by a second order Runge-Kutta algorithm. 
We consider the system size L — 8 and the self- 
inductance C ^ I. Depending on values of Iq and lu 
the number of samples used for the disorder-averaging 
ranges between 15 and 40. The number of integration 
steps is chosen to be 10^ — 5 x 10^. 
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FIG. 1. a) Upper panel: the temperature dependence of 
V^/Io for L = 8, £ = 1 and uj = 0.001. The open triangles, 
squares and hexagons correspond to Jo ~ 0.1,0.05,0.02 and 
0.01. The arrows correspond to Tp = 0.8 and Teg = 0.286, 
respectively, b) Lower panel: the size dependence of Kj//o 
for Jo = 0.05, £ = 1 and lo — 0.001. The number of time steps 
is equal to 10^. The results are averaged over 15 - 40 samples. 

The temperature dependence of the linear resistivity 
po = ^i!j/Io for different values of /q is shown in Fig. 
1 (upper panel). At low temperatures we observe a 
weak dependence on Iq, but for currents small enough 
Po becomes independent of current. From the lower 
panel of Figure 1 it is clear that the V^/Iq becomes size- 
independent for L > 10. Thus, the linear resistivity is 
nonzero for all temperatures T > in the thermody- 
namic limit. This is in good agreement with the result 
that 1jJjie_LJ(l) symmetry is not broken in the chiral glass 
statejErtil and therefore there is no superconductivity for 
any finite T. We note that a similar result was obtained 



for the vortex glasSpState when the magnetic screening is 
taken into account.EEl 

In Fig. 2 we analyze the non-linear resistivity p2 — 
—4:V^^{T)/Iq. We find that it has a negative maximum 
at a temperature Tp. This characteristic maximum de- 
pends on Iq, but we can fit its position in temperature 
at Tp = 0.8 ± 0.05 for all values of Iq presented in Fig. 
2. The arrow in Fig. 1 also indicates the position of the 
temperature Tp. We see that for T ^ Tp the linear re- 
sistivity pq is independent of current for a large range 
of currents Iq. On the other hand, below Tp the resis- 
tivity Pq is current dependent for an intermediate range 
of Jo and only for very low currents po becomes current 
independent. 

We identify Tp to correspond to the intergrain ordering 
transition temperature above which the thermoremanent 
magnetization disappears in the experiment of Ref. \V^. 
In order to verify this, we study in Fig. 3 the magneti- 
zation at a finite magnetic field / ~ HS/(pQ = 0.1. We 
show both the zero-field cooling (ZFC) and field cool- 
ing (FC) curves.a We can see that Tp is the temperature 
where there is an onset of positive magnetization, i.e. the 
paramagnetic Meissner effect starts to be observed. On 
the other hand, the irreversibility point occurs at tem- 
peratures lower than Tp, and its position is dependent on 
the heating or cooling rate. It should also be noted that 
above Tp the real part of the linear magnetic susceptibil- 
ity vanishes (see Fig. 18 from Ref. |^). 
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FIG. 2. The same as in Fig. 1 but for \Vi^{Tp)/li\. 

The results presented in Fig. 1, 2 aad 3 are in good 
agreement with the experimental data.EZI From this-point 
of view our findings and the experimental repiltslHI may 
seem compatible with the chiral-glass picture.EiJ However, 
Tp is remarkably higher than the chiral glass temperature 
Ten obtained previously (for C = l,Tcg = 0.286, see Ref. 
no) . Then we conclude that the peak of p2 has no rela- 
tion to the chiral glass phase transition. Thus, Tp just 



separates the normal state phase from a "chiral param- 
agnet" where there are local chiral magnetic moments. 
These local moments can be polarized under an external 
magnetic field, an therefore one can observe the param- 
agnetic Meissner effect under a low external field below 
Tp. At a lower temperature, collective phenomena due to 
the interactions among the chiral moments will start to 
be important, leading to the transition to the chiral glass 
state. This last transition should show in the nonlini 
chiral glass susceptibility wich should diverge at Tcg,Ej 
The chiral glass transition may also be reflected in the 
irreversibility point in the FC and ZFC magnetizations. 
Although our model is different from the corresponding 
gauge glass modeled one can expect that here the screen- 
ing spoils any glassy phase except the chiral glass. The 
linear resistivity is, therefore, nonzero for finite temper- 
atures. 
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we have calculated the specific heat, C«, which is de- 
fined as porportional to the energy fluctuations, C^ = 
{{dE)'^)/kBT^- The results are shown in Fig. 5. There is 
a broad peak in Cy located at Tp and well above Teg. Sim- 
ilar to the spin glass case where the peak of specific heat 
is positioned-higher than the critical temperature to the 
glass phase ,c3 we conclude that Tp does not correspond 
to a phase transition to a long-range ordered phase. 
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FIG. 4. The temperature dependence of V3^{T)/Io. 
The solid triangles, squares and hexagons correspond to 
iu = 0.001,0.0005 and 0.0002, respectively. L = 8, £ = 1 
and Jo = 0.1. The results are averaged over 15 samples. 



FIG. 3. The temperature dependence of the magnetization 
m in FC and ZFC regimes for the d-wave superconductors. 
L = 8 and £ = 1. The results are averaged over 25 samples. 

Our calculation of the nonlinear ac resistivity p2 is a 
non-equilibrium calculation at a finite frequency uj and 
finite ac current amplitude Iq. Therefore, one should be 
concerned on the finite lo and finite /q effects. In partic- 
ular, one may ask if it is possible that the tempereture 
Tp of the peak in p2 will tend to Teg in the limit w — > 0, 
Iq ^ 0. We have carefully studied this possibility. Fig. 4 
shows the temperature dependence of p2 for various val- 
ues of UJ and Iq = 0.1. From Fig. 2 and 4 it is clear that 
the position of Tp depends on Iq and w very weakly. It 
is, therefore, unlikely that Tp tends to Teg as a; — » and 

In accordance with the experiments of Yamao et aL,E_l 
the negative maximum of V^^{Tp)/ 1^ is shown up. Fur- 
thermore, the height of peaks of \V^^{Tp)/P\ increases 
with the decrease of lo and saturates at small frequen- 
cies (see Fig. j^). Such tendency was also observed 
exp er iment allyO 

In order to get more insight on the nature of Tp 
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FIG. 5. The temperature dependence of Cv obtained by 
Monte Carlo simulations for L = 8 and C = 1. The results 
are averaged over 20 samples. The error bars are smaller than 
symbol sizes. 

A more convincing conclusion about the nature of the 



peak in the nonlinear susceptibility should be obtained 
from a finite size analysis. Fig. 6 shows the dependence of 
max\V^^/lQ\ on the system size L for /q — 0.05 and cu = 
0.001. Clearly, the height of the peak docs not diverge 
as L — + oo. In other words the peak in the nonlinear 
resistivity does not corresponds to a phase transition in 
the thermodynamic limit. 




FIG. 6. The dependence of the maximal values of 
\VLiTp)/Io\ on the system size L. Iq = 0.005, £ = 1 and 
uj = 0.001. The results are averaged over 15 - 40 samples. 



(0 













1 


cj = 0.001 




1 


- 









- 


1 . ^, ,^ 


- 


1 

o 


1 


1 1 


• 



a = 1.1 ± 0.6 giving more weight to small values of Iq. 
So within the error bars our esti mate of a agrees with 
that obtained by the experiments. tZHl3 



IV. DISCUSSION 

In Ref. |lj it was argued that the peak of the nonlinear 
resistivity was a signal of the transition to the chiral glass 
state. The value of Tp obtained in our simulations is, 
however, considerably higher than the chiral glass tran- 
sition temperature. Teg. We conclude that the peak of 
P2 is not related to the transition to the chiral glass. 
Tp is found to coincide with the point for the onset of 
the paramagnetic Meissner effect, where the magnetiza- 
tion becomes positive. In thisjjespect, our result agrees 
with the experimental result .EZl We interpret Tp as the 
crossover temperature from the normal state phase to a 
"chiral paramagnet" in which there are local chiral mag- 
netic moments induced by the 7r-junctions. As the tem- 
perature is lowered down the system would have a phase 
transition from the chiral paramagnetic phase to the chi- 
ral glass state. At this critical point p2 does not show 
any particular feature. Furthermore, we found that the 
linear resistivity is always finite at T > due to screen- 
ing effects, and therefore there is no superconductivity in 
the random tt junction model. 

Xu conclusion, the experimental results of Yamao et 
alxli can be reproduced by the XY-like model for the 
d-wave superconductors. Contrary to the speculation of 
Ref. n^ we expect that Tp does not correspond to the 
chiral glass transition. 
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FIG. 7. The dependence of the maximal values of 
\VLiTp)/Io\ on Jo. L = 8,£ = 1 and uj = 0.001. The re- 
sults are averaged over 15 - 40 samples. 

Fig. 7 shows the log-log plot for the dependence of the 
maximal values of | V^^ {Tp )/I^\ on /q for a fixed frequency 
Lj = 0.001. One can fit max IK/, ,(T„)//^| - //? with 
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